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ABSTRACT 


An approximate analytical solution is obtained for 
the temperature distribution in an absorbing gas layer 
bounded by two flat surfaces radiating at different tem- 
peratures. The corresponding expression for the net 
radiation exchange between the two surfaces is derived 
and numerical solutions obtained for gas layers up to 10 
optical thicknesses in depth. Similar results are derived 
for the special case of a gas layer at uniform temperature, 
and the net radiation exchange in this system is found to 
agree well with that in the previous for gas layers less 
than two optical thicknesses in depth. 


I. INTRODUCTION 


A problem in radiative transfer of increasing importance to the calculation of the thermal balance in 
hi gh- temperature propulsion devices is the net radiation exchange between two surfaces at different tem- 
peratures separated by an absorbing gas. Various specialized aspects of this general problem have received 
the attention of investigators from time to time during the past century. Some of the methods used in these 
analyses have been summarized by Jakob (1). Of particular interest are the analytical methods of Wohlenberg 
(2) and Christianson (3) for obtaining the net radiative transfer in systems wherein the absorbing gas is at a 
uniform temperature. More recently, the detailed problem of determining the steady-state temperature distribu- 
tion in the gas, as well as the net energy transfer, has been treated by Usiskin and Sparrow (4) who obtained 
numerical integrations of the appropriate integral equations. 


1 


JPL Technical Report No. 32-197 


The present treatment examines some of these same aspects of the problem by means of a dual 
approach. In the first portion of the analysis, the methods of Wohlenberg and Christianson are combined so 
as to handle the case of two grey surfaces separated by an absorbing gas at uniform temperature. In the 
second portion, the problem considered by Usiskin and Sparrow is reviewed, and an approximate analytical 
solution is derived which yields closed-form expressions for the temperature distribution in the gas and the 
net radiative transfer across the system in terms of the relevant physical parameters. 
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II. ABSORBING GAS AT UNIFORM TEMPERATURE 


The problem considered here is that of determining the radiation exchange between two flat grey 
surfaces separated by an absorbing gas at uniform temperature. The absolute temperature of one surface is 
Tj and of the other, T^, where T ^ > T ^ It Is further specified that the surface at temperature T- has 
emissivity and that this is equal to its absorptivity. The gas has absorptivity a. to the radiation emitted 
from the surface at temperature T-. The radiation balance relations are stated in terms of three fundamental 
quantities: 


V. s cr e. T f = grey-body radiation (energy/area) from surface at temperature T-; (1) 

cr is the Stephan — Boltzmann constant. 

W.j = radiation /area from medium i absorbed by medium /. (2) 

q. . = direct radiative transfer/area between media i and j. (3) 

*»/ 

A. Radiation from Gas 


The first step in the analysis is the derivation of an expression for the radiation from the gas. For 
this purpose we require such quantities as W «, the radiation emitted by the gas which is absorbed by the 

g 2 

wall at temperature T ^ The quantity W 2 consists of two parts, that portion of the gas radiation which goes 
directly to wall 2, and that which goes first to wall 1 and, after reflection, is eventually absorbed by wall 2. 
The application of Christianson 9 s method of multiple reflections yields for the former: 


W 


+ (1 - c 9 ) (1 - a 9 ) (1 - € 2 ) (1 - a, ) e, 




+ 


[(1 - e 2 )(l - a 2 )(l - e 1 )(l -a^]' 


£ 2 + 


e 2 ^ 


e 


l-V 


where 


1? = (1 - Oj) (1 - Oj) (l - SjXl- e 2 ) 


( 4 ) 
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is the fraction of the radiation leaving either solid surface which is reflected once by each; the last equality 
in the above equation follows by virtue of the fact that V< 1. By applying this approach also to the 
radiation component which goes first to wall 1, is reflected and subsequently absorbed by wall 2, one obtains 

e i )(1 " a i> 

(1 - l) 


Adding these expressions yields W 


gr 


W 


s 2 


[1 + (1- ffjMl-aj)] 

1-V 


(5) 


A related quantity of interest is the radiation emitted by wall 2 which is ultimately absorbed 

by the gas. By a similar line of reasoning to that applied above, it is easily shown that 


W„ 


2g 


1 - Tj 


[ a 2 + Oj (1 - £j) (1 - a 2 )] 


( 6 ) 


It is of interest to mention that the radiation emitted by wall 2 and subsequently reabsorbed, IF 22 , * s 8* ven by 


W 

"22 ” 


lT 2 (l-a l )(l-a 2 )(l- £ x ) e 2 

1-V 


(7) 


and, of course. 


^2 ~ ^22 + ^21 + ^ 2 g 


( 8 ) 


Corresponding expressions for W ^ and W ^ are obtained by simply exchanging indices in the above 


relations. 
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An explicit solution for V (which yields the equilibrium temperature of the gas in this system) is 
obtained by requiring that the net radiation passing from wall 1 to gas is equal to that from gas to wall 2. 
Using the notation of definition (3), the latter is given by 


R 


g, 2 


e 2 W [l + d-CjHl - a j)]-r 2 [cu + aM-e )(i-a,)] 

w _ w = 5 

W g 2 W 2g 

I~V 


(9) 


and the former by —q , which may be obtained immediately from (9) 

g » 1 * 

steady state. 


by exchanging indices. Now, at 


= q 


g. 2 


( 10 ) 


from which follows 


V 2 [“2 + “i (1 - e x)(l - “ 2 )] +1^1 [ 0 ^ + a 2 (l - e 2 Hl - 
e 2 [l + (1- fijHl-Oj)] + 6j [l + (1- e 2 )(l-a 2 )] 


(ID 


B. Not Radiation from Wall 1 plus Gas to Wall 2 


The combined net radiation from wall 1 and gas to wall 2 is given by 


? l+g.2 


W 12 + W g2 


(W, 


21 


W 2g ) 


( 12 ) 
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where W ^ and W are obtained from Eqs. (5) and (6), and 


W 


12 


V 1 d-V e 2 


1 - 7 ] 


w 


21 


*^2 ^ ~ a 2 ^ € 1 


1-17 


( 13 ) 


Substitution of these expressions into Eq. (12) yields 


<? 1+g2 = — {^(1-a,) e 2 + e 2 ir g [l + (1 - ej (1 - Oj)] 

1 - 7 ] 


- IT 2 [(1 - a 2 ) €j + a 2 + Oj (1 - €j) (1 - a 2 )] } 


(14) 


C. Direct Radiative Exchange Between Wall 1 and Wall 2 


The direct radiative exchange between 


walls 1 and 2 is obtained immediately from results (9) and 


(14); thus, 


q \,2 “ q l+g,2 ~ q g,2 ~ W \2 W 21 


(15) 


Substitution from Eq. (13) yields 


9 


1,2 


W \ e 2 a - Oq) - H 2 (1 - a 2 ) 

1 - 7 ] 


(16) 
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# 


D. Limiting Cases 

Three limiting cases for the above results are of particular practical importance. 

Consider first the case when both surfaces are black-body radiators; then 6^ = = 1. The 

corresponding expression for the radiation from the gas is obtained from Eq. (11). 


V g 2 ^2^2 + ’’V 


Hie net radiation exchange between wall 1 and gas and wall 2 then reduces to 


(17) 


q l + g,2 ~ ^ °T^ V l + V g~ W 2 


(18) 


which is in agreement with the result given by Jakob! The substitution of the form (17) into (18) yields 


q l+g,2 - ft — °i) *1 - d - — “ 2 ) V 2 


(19) 


Finally, one can also obtain the appropriate result for the net exchange between walls 1 and 2 from (16). 


7l,2 “ d “l) ^1 ~ (1 ~ a 2) ^2 


( 20 ) 


* Ref. 1, Vol. II, p. 132. 
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The second limiting case of interest is that of an entirely opaque gas, that is, = a 2 = 1* three 
principal results in this case reduce to 


\ = 


V 2 +W l 


e i +e 2 


( 21 ) 


^1 + g »2 = £ 2 % 


= e n W_-W 2 = 




£ 1 + £ 2 


( 22 ) 


where the last equality is obtained by the application of Eq. (1) and (21). The direct exchange between walls 
1 and 2 is then 


? 1,2 = 0 


(23) 


as expected, since the gas layer completely shields one wall from the other. 

The last special case is that of a completely transparent gas, that is, = 0. Then, 


- 0 


(24) 


^l+g,2 " 


cr e i « 2 (r}-r{) 

1 - (1 - €j) (1 - e 2 ) 


(25) 


again using Eq. (1). This last expression is the familiar form reported in the literature (see Refs. 5 and 6). 
The corresponding expression for q j 2 is identical to Eq. (25), as is to be expected, since the gas layer 
cannot influence the energy balance by virtue of its transparency. 
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III. ABSORBING GAS WITH EQUILIBRIUM TEMPERATURE PROFILE 


The simplified physical problem treated in the preceding section is now generalized to account for 
the presence of the gas layer as an emitting and absorbing medium whose internal temperature at any local 
station is determined by the net energy exchange with all other radiating elements in the system. In carrying 
out the analysis we introduce the following simplifying assumptions: (1) the two solid boundaries radiate 
as black bodies; (2) the gas radiates as a grey body; and (3) heat transfer by conduction is negligible. 

An immediate consequence of assumption (2) is that the absorption characteristics of the gas layer 
may be represented by a single coefficient of absorption (i.e., the logarithmic decrement of radiation) which 
is independent of wavelength. For further simplification, we assume that this coefficient is also independent 
of temperature. By virtue of assumption (3), the energy balance in the system may be specified entirely in 
terms of the radiative balance. The appropriate equation is obtained from an analytical statement of the 
energy balance in a differential volume element in the gas and yields a linear integral equation for the 
temperature distribution within the gas. This equation has been derived by Usiskin and Sparrow (4) and is 
given as Eq. (15) in the reference paper. 

A. Temperature Distribution 

The purpose of the present section is to develop an approximate analytical solution for the 
“temperature” (actually radiative power) distribution which involves directly the essential physical para- 
meters of the system, and subsequently to apply this result to the calculation of the net radiative transfer 
across the gas gap. The governing energy balance equation of Usiskin and Sparrow is taken over directly 
for this analysis, and it may be shown, upon introducing the generalized exponential integrals and some of 
their properties, that this equation may be written in the following form 


2 W g (x) 


Six) 

2 k 


+ ^2^2 + W l E 2 (k(1 - x)) + k f W g is)E l (k(x - s)) ds + k /V {s)E l (k is - x)) ds 


( 26 ) 


where 0 <_ x < 1 is the fractional distance across the gas layer measured from the cold (Tg) wall. The 
quantity 4k ix)dx represents the energy radiated per unit time from a gas layer of thickness dx and unit 
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cross-sectional area (2). A special case of W (x) is the familiar expression cr T 4 (x), where T (*) is the 

6 & & 

temperature of the gas. The symbol W - is defined as in the preceding analysis. The function S (%) represents 
the internal energy sources in the gas (with units of power/volume); the parameter k is the optical thickness 
of the gas layer and is defined k = kL; k is the absorption coefficient of the gas and L is its thickness. The 
E n (x) are the generalized exponential integrals; these are defined (7), 


E n (*) ^ 


„n-l 


r 


ds 


Equation (26) states simply that the power radiated from a unit volume of gas at x (the term on the 
left-hand side) is equal, at steady state, to the radiated power it absorbs from all other elements. The first 
term on the right-hand side accounts for the internal sources; the second term gives the portion of the 
radiation emitted by the cold wall which is absorbed by the gas volume, and the third gives that portion from 
the hot wall. The third and fourth terms account for radiation from the rest of the gas layer which is absorbed 
in the unit volume about x. 

This energy balance equation may be considerably simplified by introducing the dimensionless 
function 0 (x): 


w (x)~ a T* 

6(x) = (27) 

<r(T\- T\) 


The substitution of this function into Eq. (26) yields 


2 6{x) 


Six) j 

+ EAk (1 - *)) + K f d(s) E, (k I x - s| ) ds 

2 kaiT\-T\) 0 


(28) 


where we have used the result 
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In this analysis attention is focussed on the solution to Eq. (28) for the case of no internal heat 
sources in the gas; thus, SW = 0. The resulting equation is in the form of the inhomogeneous Fredholm 
equation of the second kind, and formal solutions can be generated by, for example, a Liouville-Neumann 
series expansion (8). An examination of the numerical results obtained by Usiskin and Sparrow reveals, 
however, that the “temperature ' 9 function 6{x) has nearly a linear dependency on x , especially for large k- 
Therefore, in constructing our approximate analytical solution we introduce the form 6{x;k) - o(k)x + c (k) 
into the integral as a zeroth-order approximation and calculate the resulting expression for Q(x; k) at the 
left. If the assumed linear form is indeed a reasonable first estimate for 0(x), then the resulting expression 
from the left should have the form 


0 (. x ; k) = a (k) x + c (k) + R (x; k) 


(29) 


where the correction term R (x; k) will be small. 

For computational purposes, it is convenient to write the source-free equation in the form 


2 6{x) = E 2 (k(1-x)) + k f 6(s) E^(k(x - s)) ds + k f 0(s) E ^(k{s - x)) ds 


(30) 


The application of the procedure mentioned above is straightforward, and leads rather easily to the form (29). 
In carrying out this calculation it is necessary to impose also two conditions on 6(x) so that the specification 
of the quantities a (k) and c (k) is complete. These are 


0 



1 

2 



(31) 


The first condition is a consequence of the antisymmetry property of the function 0 (see Appendix for proof). 
The second condition implies that R 1 (1/2; k) - 0. The substitution of the linear form for 6 into the right-hand 
side of Eq. (30) yields the result: 


o(k) 

e ~ KX _ e -K(l-%) 


1 

c{k) H a(/c)x 

E 2 (kx) + 

1 

1 - c(k) a(x) (1 + x) 

£ 2 (k(1 - *)) 

2k 

1 

- 


2 


2 



( 32 ) 
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4 


In performing the calculation the following property of the E n (*) is required: 


J sE l {s)ds = xE 2 (x)-E 3 {x) 

0 2 


The application of relations (31) yields 



c{k) = — [ 1 - a (k) ] 

2 


(33) 


The resulting expression for 6(x;k) given in Eq. (29) (see Fig. 1 ) when compared graphically with the 
numerical solution obtained by Usiskin and Sparrow is found to be accurate to less than a percent. 

The comparison is especially good for large k; when k > > 1, the expressions for a and c have the asymptotic 
forms: 


2 

a (k) ~ 1 - — 

K 



(34) 


i 

It may be observed from the form of the function (29) that further iteration will yield an integral 

equation for the correction term R(x;k). Because of the antisymmetry property of 0 , the next-higher-order 

I 

approximation may be constructed by introducing an x ^ term into (29). A continuation of this procedure 
will generate a power series in odd x which will converge, thereby providing an increasingly accurate 
analytical expression for 6. 

B. Net Radiative Transfer 

Of practical interest is the calculation of the net radiative transfer between the hot wall and gas 
and the cold wall. This corresponds to the quantity < 7 , 9 determined in the previous section for the case 
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of an isothermal gas [see Eq. (14)] . Again, the appropriate expression has been derived by Usiskin and 
Sparrow. 



FRACTIONAL DISTANCE FROM COLD WALL, x 


Fig. 1. Approximate solution for temperature distribution function 


By introducing the E n (x) functions and combining various terms, their result may be written in the form 


Q{k) — 2E«{k) + 2k J 0(s) E tyitts) ds 
6 0 z 


(35) 


where 


qU) 

QM = 

cr(T\- T*) 


(36) 


and q (k) is the net energy transferred per unit area and time to die cold wall. It is easily shown that Q(k) 
has the following properties: 


lira Q(k) - 1 
k - 0 


lira Q{k) = 0 (37) 

k-0 


13 



JPL Technical Report No. 32-197 


as is to be expected on physical grounds. The function Q (k) has been computed by Usiskin and Sparrow for 
the interval (0 <_ k < 2) using the numerical solution for 0(x) mentioned previously. Of particular interest 
to this discussion is the comparison of their result with that obtained using the analytical form for 6 
suggested in the preceding section. If one neglects entirely the correction term /?, then direct substitution 
into Eq. (35) yields 


Q(k) 




a ( k ) 


2 a (k) e 


3 K 


aW 


1 - 


f|W 


(38) 


It may be shown that this result also satisfies the limits (37). A graphical comparison with the exact 
numerical solution is shown in Fig. 2. Evidently the use of the linear form for 0 yields an estimate for Q 
which agrees to within about 10 percent in the interval (0.1 < k < 1); even better agreement is obtained 
outside the interval. 

A final comparison of some practical value may be drawn from the result obtained for the net 
radiation exchange Q{k) using the relations derived for the uniform temperature gas system discussed in the 
preceding section. The appropriate expression is given in Eq. (19), which for the present system may be 
written 


QU) 


? l + g,2 


W - w 

w l w 2 


- 1 


1 

— a (k) 
2 


where - a. From the definitions of a and k, it follows that 


a (k) = 1 - e ~ K 


(39) 


Thus, for the uniform temperature gas, 


<?(«) = — (1 + e" K ) 
2 


( 40 ) 
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This function is shown in Fig- 2, and the agreement with the variable temperature case is reasonably good 
up to two optical thicknesses of gas layer. 



Fig. 2. Net radiative transfer to cold wall 

It is interesting to note the difference in the limits for the two cases when k -> «>• The result for 
the case with non-uniform gas temperature was discussed previously [see Eq. (37)] . The result for the 
uniform temperature gas case is also explainable on physical grounds. In the limit of a very opaque gas, 
it is to be expected that the net radiative transfer to the cold wall will be simply the difference between 
the radiation emitted by the gas and the cold wall; thus, one can write immediately 


lim Q (k) 
00 


W g (oo)-W 2 



(41) 


where W g (°°) denotes the radiation emitted by an opaque gas (i.e., a-> 1, k -* <»)- If the expression (21) is 
used for this quantity, along with = 1, then it follows that Q (/<)-► 1/2 as k -> 
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APPENDIX. The Antisymmetry Property of the Function 9(x). 


Consider the function d(x) which satisfies the integral equation 


[cf. Eq. (30)] 


2${x) = £„(k(1 - x)) + k f 1 9(s) E, (k|* - s | ) ds 

i 0 1 


0 < * < 1 (A-l) 


Transform to the function <f>(x), where 


e (*) = e{x) + 


(A-2) 


Substitution into Eq. (A-l) yields 


2 <£(*) + 1 = £<,(k(1 - x)) + k f 1 <?b (s) E, U|* - s| ) ds + — E.(k\x - s\) ds (A-3) 
z 0 2 0 

also 

1 K 1 

2 d) (1 - *) + 1 = E<,(kx) + k f <i>(l - s) E , (k \ x - s\) ds + — f E ,(k\x - s\) ds (A-4) 
* 0 2 0 


If a new function f{x) h <£> (x) + cj> (1 - x) is introduced, then the addition of Eq. (A-3) and (A-4) reveals that 
f(x) satisfies the integral equation 


2 f(x) = K f l f(s) E,(k\x - s\) ds 
0 1 


(A-5) 


wherein use has been made of the relation 


K 




E, («|* 


s\) ds 


2 - E 2 (kx) - E 2 (k{ 1 - x)) 


(A-6) 


* 


The author is indebted to C. Solloway for assistance in demonstration of this proof. 
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Clearly, f(x) - 0 is a solution to Eq. (A-5). This yields the property 

0(*) + 0(1-*) = 0 (A-7) 

which is in fact the statement of antisymmetry for the function 0(*); i.e., 0 (*) is an antisymmetric function 
of the argument [* — (1/2)] - It is obvious then that 6(x) is antisymmetric about the point (1/2, 1/2) in the 
0, * -plane. 

To prove that f{x) ~ 0 is a unique solution to Eq'. (A-5), suppose that there is some non-zero 
solution /(*). Then, 


2\f(x)\ < \ f\ K J 1 EAK\x-s\)ds 
max 0 

<r\f\ 

max 

where r < 2, which follows from Eq. (A-6). Furthermore, 

max 

where £ s r/2 and 0 < £ < 1. Since this inequality holds for all *, it follows that 

l/l <£\f\ <1/1 (A-8) 

max max max 

which is a contradiction. Therefore f(x) = 0 is a unique solution, and 0(*), and likewise 6{x), are 
antisymmetric functions. 


17 


JPL Technical Report No. 32-197 


REFERENCES 


1. Jakob, M., Heat Transfer , John Wiley and Sons, Inc., New York, Vol. I, 1949, Vol. II, 1957. 

2. Wohlenberg, W. J., Eng. Bull. Purdue University, Vol. 24, No. 4a, Engineering Experiment Sta. Res. 

Ser. No. 75, August 1940. 

3. Christianson, C., Wiedemanns Ann. 19, 267, p. 1883. 

4. Usiskin, C. M. and E. M. Sparrow, Int. }. Heat Mass Transfer, 1 , Pergamon Press, New York, 1960, 
pp. 28-36. 

5. Jakob, M. and G. A. Hawkins, Elements of Heat Transfer and Insulation, 2nd Ed., John Wiley and Sons, 
New York, 1950. 

6. Condon, E. U. and H. Odishaw, Editors, Handbook of Physics, McGraw-Hill Book Co., Inc., New York, 
1958, p. 5-76. 

7. Case, K. M., F. de Hoffman and G. Placzek, Introduction to the Theory of Neutron Transport , Vol. I, 
Los Alamos Scientific Laboratory, New Mexico, 1953 (available through U. S. Government Printing 
Office, Washington, D. C.) 

8. Morse, P. M., and H. Feshbach, Methods of Theoretical Physics, Vol. I, McGraw-Hill Book Co., Inc., 
New York, 1953, Ch. 8. 


18 



